A Prediction of the B* mass in full lattice QCD 
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By using the Highly Improved Staggered Quark formalism to handle charm, strange and light 
valence quarks in full lattice QCD, and NRQCD to handle bottom valence quarks we are able 
to determine accurately ratios of the B meson vector-pseudoscalar mass splittings, in particular, 
(m(B*) — m(B c ))/(m(B*) — m(B s )). We find this ratio to be 1.15(15), showing the 'light' quark 
mass dependence of this splitting to be very small. Hence we predict m(B2) = 6.330(7)(2)(6) GeV 
where the first two errors are from the lattice calculation and the third from existing experiment. 
This is the most accurate prediction of a gold-plated hadron mass from lattice QCD to date. 



INTRODUCTION 

Particle physicists now have long familiarity with the 
low-lying spectrum of bb (Upsilon) and cc (psi) mesons 
but they nevertheless continue to provide a very impor- 
tant testing ground for our understanding of strong in- 
teraction physics. The similar be (B c ) system, on the 
other hand, is largely unexplored territory so predictions 
of the meson masses are very valuable. These predic- 
tions need to be as accurate as possible (and with an 
error budget) to provide stringent tests of QCD. Lattice 
QCD is clearly one of the best ways to do this, now that 
accurate calculations including the full effect of u, d, and 
s sea quarks inside hadrons are possible p]. It has al- 
ready provided successful predictions of the pseudoscalar 
rit, mass [5] (with a 14 MeV error) and the B c mass [3] 
(with a 20 MeV error), both subsequently seen by exper- 
iment. Here we give a prediction for the vector B* mass 
through the mass difference between the B* and the B c . 

Mesons composed of valence heavy (b and c) quarks 
arc relatively simple because they are nonrelativistic sys- 
tems and a potential model may be expected to work 
reasonably well (see, for example, [HOE]). This is es- 
pecially true for the T system where v\ rs 0.1 (in units 
of c 2 ). It is less true for charmonium where u 2 ps 0.3 
and so relativistic corrections are much larger there. The 
ground state hyperfme (vector-pseudoscalar) mass split- 
ting is such a correction, but is given in leading order 
perturbation theory by a simple formula since the S.S 
potential is proportional to 6 3 (f). 

AM= 32^(0)1* (1) 
9mim2 

where mi and mi are the masses of the quark and anti- 
quark and ip(0) is the wavefunction at the origin from the 
potential model. Using this formula to calculate the split- 
ting will have a systematic error at 0(v 2 ) i.e. 30% in cc, 



10% in bb and 20% in be. However, the cc hyperfme split- 
ting has been used in the past to fix the effective value of 
a s in eq.[T|and then that 30% error affects all subsequent 
calculations. A larger problem, perhaps, is the variation 
in results between different potential models tuned to the 
spin-independent spectrum. This is because that spec- 
trum does not in practice constrain the wavefunction at 
the origin at all strongly. The mass splitting between B* 
and B c can vary in the range 40-90 MeV |U \5\ |B] be- 
tween different potentials, which makes it hard to decide 
a 'central value' and error budget. 

The reduced mass in the B c system is roughly one half 
that of bb and 1.5 times that of the cc. Then w 2 « 0.05 
in B c but v 2 ~ 0.4 — 0.5, which makes a nonrelativistic 
treatment worse in principle than for charmonium. An 
alternative approach is to treat the B c as a 'heavy-light' 
system using ideas from HQET but, for example, it is 
difficult to estimate the light quark mass dependence of 
the 1/toq operator giving rise to the hyperfme splitting, 
limiting again the accuracy in the prediction. 

Lattice QCD, on the other hand, can provide very 
stringent tests of QCD from the hadron spectrum, in 
which all sources of systematic error can now be tested 
and quantified [1 . The only parameters are those of 
QCD itself (a quark mass for every flavor and a coupling 
constant) and impressively accurate results in agreement 
with experiment can be produced for the whole range of 
gold-plated hadron masses known experimentally. Our 
previous prediction of the B c mass [3] dates from the rel- 
atively early days of full lattice QCD calculations and is 
now being improved. We have since developed a much 
more accurate method for handling charm quarks within 
lattice QCD and that has enabled a determination of the 
B c mass with smaller systematic errors [7] . This method 
also allows an accurate prediction of the B* and we de- 
scribe that calculation here. 
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TABLE I: Ensembles (sets) of MILC configurations used with 
gauge coupling f3, size L 3 x T and sea masses (x tadpole 
parameter, no) rn^ q and m^ 9 . The lattice spacing values 
in units of ri after 'smoothing' are given in column 3 [8]. 
Column 8 gives the number of configurations and time sources 
per configuration that we used for calculating correlators. On 
set 5 only half the number were used for light quarks. 
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TABLE II: Parameters for the valence quarks. aM® is the b 
quark mass in NRQCD, and uol is the tadpole-improvement 
parameter used there [5] . We use stability parameter [2] n = 4 
everywhere. Since NRQCD quarks propagate in one direction 
in time only we improve statistics by generating propagators 
both forwards in time (for T/2 time units) and backwards in 
time from each source. Columns 4, 6 and 7 give the charm, 
strange and light bare quark masses for the HISQ action. 1 + e 
is the coefficient of the Naik term in the charm case [12] . 



LATTICE QCD CALCULATION 

From above it is clear that an optimal lattice QCD ap- 
proach to the B c is to combine a nonrelativistic method 
for the b quark with a relativistic one for c. Here we 
use Lattice NRQCD for the b, developed over many 
years [5J [TO] E] to provide accurate bottomonium spec- 
troscopy [2 by including spin-independent terms through 
0(v^) and leading spin-dependent terms with discretisa- 
tion corrections through 0(a 2 ). For the c quark we use 
Highly Improved Staggered Quarks (HISQ) [15], a fully 
relativistic discretisation of the Dirac action which is ac- 
curate enough to handle c quarks because it is fully im- 
proved through 0(a 2 ) and also has the leading (m c a) 4 
errors removed. This enables us to use the same lattice 
QCD action for charm, strange and light quarks (we take 
m u = rrid). This approach, as we shall see, enables us 
to cancel some systematic errors between the B c system 
and the B s system and obtain the hyperfine splitting in 
the B c as a multiple of the experimentally known split- 
ting [13] in the B s system. 

We work with five ensembles of gluon field configura- 
tions provided by the MILC collaboration. These include 
the full effect of u, d and s sea quarks using the improved 
staggered (asqtad) formalism and are available with large 
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FIG. 1: The ratio R c of the B c and B s hyperfine splittings 
(eq. [3| as a function of lattice spacing, a, from full lattice 
QCD. Our continuum extrapolation is also given and the re- 
sult at a = 0. The lighter points and line give the equivalent 
points for 71; along with the experimental value [13] . 

spatial volumes (> (2.4fm) 3 ) and at multiple values of 
the light sea masses (using m u — m^) for a large range 
of lattice spacing values. We use configurations at three 
values of a between 0.15 fm and 0.09 fm with parameters 
as listed in Table [I] On each configuration in the ensem- 
ble we generate b quark propagators using NRQCD and 
c, s and I quark propagators using HISQ. The parame- 
ters of the valence quarks are given in Table |TT] The b 
quark mass is tuned to give the correct T mass |14] and 
the charm, strange and light masses are taken from [15] . 

The b quark is then combined in turn with each of the 
other three with appropriate spin matrices to make pseu- 
doscalar or vector mesons. To increase statistics we gen- 
erate propagators from sources at several different times- 
lices per configuration (see Table [T]) . We also use a ran- 
dom wall source for the quarks [15], taken as a set of U(l) 
random numbers at each point on the source time slice. 
This mimics multiple sources across a timeslice when the 
propagators are paired up, improving statistics further. 
For the NRQCD propagators, as well as a local source, 
we also need 'smeared' sources [5] chosen to improve the 
overlap with the ground state in the meson correlator. 
Exponentially growing noise is a problem in the B system 
(particularly as the lighter quark mass becomes small) 
and smearing enables us to extract an accurate ground 
state energy from the correlator at smaller time separa- 
tion from the source than otherwise [TB] • We use a Gaus- 
sian form for the smearing function with radius 2a and 
4a. These various sources for the NRQCD quark must 
be combined with the random wall described above. In 
addition the NRQCD quark source must now include the 
matrix that converts spinless staggered quarks into naive 
quarks for combination with 2-spin NRQCD quarks in 
an adaption [7] [16] of the standard method of combining 
heavy quarks with staggered quarks [17] . 

We fit our 3x3 matrix jTH] of B meson correlators 
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FIG. 2: The spectrum of 'gold-plated' mesons from HPQCD calculations. Results are divided into those used to fix the 
parameters of QCD (4 quark masses and a coupling constant); those which are postdictions [21 1121 FT5] and those, like the B* 
described here, which are predictions J2j [3] . 



using the standard Bayesian method [H5] to a sum of 
exponentials, including oscillating parity partner states 
as: 



-Ek(t—to) 



(2) 



Ncxp-l 



(t-*o) 



fe'=0 



where i,j index different smearing radii. We look for 
stability in the fits and their errors as a function of iV exp 
for ground state energies, Eq. 

Eq is not the meson mass but contains an energy 
shift due to the non-relativistic treatment of the b 
quarks [HI EH E] • The shift cancels in the mass differ- 
ence between states with the same NRQCD quark con- 
tent. Thus the B*—B q splitting is obtained directly from 
A q = E (B*) — E (B q ). Because errors are strongly cor- 
related between similar quantities calculated on the same 
ensembles we fit B q and B* correlators simultaneously to 
the form above and determine A q directly from the fit. 



Values arc given in Table III for q — I, s, c. 

The terms from the HISQ action that contribute to 
the hyperfine splitting are hidden inside the discretisa- 
tion of the Dirac covariant derivative. Because HISQ is a 
relativistic action, these terms will automatically be cor- 
rect in the a — > limit. The spin-dependent term in the 
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1 0.0318(78) 0.0311(37) 0.0324(2) 1.02(27) 1.04(12) 

2 0.0374(35) 0.0359(21) 0.0326(3) 1.04(11) 0.908(53) 

3 0.0306(54) 0.0287(19) 0.0268(2) 1.06(19) 0.934(62) 

4 0.0245(68) 0.0261(27) 0.0271(4) 0.94(26) 1.04(11) 

5 0.0177(35) 0.0190(14) 0.0220(6) 0.93(19) 1.158(91) 







1.00(23) 1.15(15) 



TABLE III: Results for the mass differences between vector 
and pseudoscalar B mesons for different light quark content 
on different MILC ensembles. A q = E (B*) — E (B q ). R q is 
the ratio A q /A a and the result of extrapolating R c and to 
a — is also given. 



NRQCD action that gives rise to the hyperfine splitting 
can instead be explicitly pinpointed as the a-B term [10] . 
This term has the correct tree-level coefficient to match 
full QCD at 0(vf) but radiative corrections beyond this 
have not been included. Hence the normalisation of this 
term, and the normalisation of the hyperfine splitting, 
have an uncertainty of 0(a s ) (« 20%). This uncertainty 
is part of the NRQCD action and hence the same un- 
certainty appears regardless of which light quark is com- 
bined with the b quark and cancels in ratios of hyperfine 
splittings. In Table |TH| we also give values for 
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and the corresponding quantity, Ri, for u/d quarks. On 
sets 1-4 Ri is given directly by a joint fit to B s and Bi 
correlators. All Ri values agree with 1 within 30% errors. 

Figure [T] shows R c as a function of lattice spacing. 
There is little dependence on the light quark mass, since 
neither the B c nor the B s contain valence light quarks 
and we do not expect strong sensitivity to the sea con- 
tent. Lattice spacing dependence is mild — the dashed 
line shows an extrapolation to the continuum limit at 
a = that can be compared to experiment. The extrap- 



olation includes 



and 



terms and allows for linear 



dependence on sea quark masses. In that limit we find 
R = 1.15(15). This, along with the results for Ri show, 
somewhat surprisingly, that the hyperfine splitting varies 
hardly at all with the mass of the lighter quark in the B 
system, up to and including charm. 

The result is backed up by the existing experimental 
results on heavy- light and heavy-strange mesons. In the 
D system the hyperfine splittings differ by only 2% be- 
tween the D s and the D^. Some of this difference may 
in fact be a result of coupled-channel effects since the 
is just above threshold for the decay to Dir, whereas the 
D* has only the OZI-disfavoured decay mode D s ir avail- 
able. The experimental situation is less clear in the B 
sector since some experimental results favour a B* — B s 
splitting very close to the B* — B (not yet differenti- 
ated into charged and neutral modes) and others favour 
a somewhat larger splitting [TS]. We use the PDG aver- 
age value of 46.1(1.5) MeV [13] for the B* - B s splitting 
because, in keeping with our result and indications from 
the D sector, this is closer to the B* — B splitting than 
the PDG fit value of 49.0(1.5) MeV [T3]. 

Our result for R gives 53(7) MeV for the B* - B c split- 
ting, where the error is statistical only. Additional sys- 
tematic errors come from relativistic corrections to the 
a ■ B term in the NRQCD action [TO] . We can estimate 
the size of these from the size of v\ in the B c (0.05) 
and the B s (=(AQCD/mb) 2 = 0.01). The cancellation 
between them leads to a 4% systematic error. Any mis- 
tuning of the b quark mass will cancel in R, and small 
mistunings of the s and c quark masses lead to a negligi- 
ble error. Electromagnetic hyperfine effects missing from 
our calculation should also be negligible (less than 1%). 



CONCLUSIONS 

Adding our value for the B* — B c splitting to the ex- 
perimental mass for the B c gives the mass of the B* as 
6.330(7) (2) (6) GeV where the first two errors are from 
the lattice QCD calculation - statistics and systematics 
respectively - and the third error is from experiment for 
the B c and the B*. The relatively small value of the 
B* — B c splitting will make it challenging to find the B* 
from its decay to B c j. 

The absence of strong dependence of the hyperfine 



splitting on the mass of the lighter quark in the B system 
is an interesting result, which has implications for other 
spin-dependent splittings in the B c system. In HQET 
language it says that matrix elements of the hyperfine 
operator are insensitive to the light quark mass, up to 
and including charm. In constituent quark model lan- 
guage, using a formula for the hyperfine splitting akin to 
that in eq. [I] the result implies that |?/j(0)| 2 varies as m q 
to cancel the m q in the denominator [4j. The amplitude 
a ioc,a from the fit in eq. [3] is proportional to ^(0) in a 
nonrelativistic approach. This does show significant de- 
pendence on the light quark mass in the B. For example, 
aioc,o(B c )/ai OC fl(B s ) « 2 [?]■ 

Finally, in Figure [2] we summarise the current status of 
the gold-plated meson spectrum as determined from lat- 
tice QCD, highlighting those meson masses which have 
been made as predictions ahead of experiment. The re- 
sult here is the most accurate prediction to date. 
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